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Magnetars are surrounded by diffuse plasma in magnetic field strengths well above the quantum electro¬ 
dynamic critical value. We derive equations of “quantum force-free electrodynamics’" for this plasma using 
effective field theory arguments. We argue that quantum effects do not modify the large scale structure of the 
magnetosphere, and in particular that the spin-down rate does not deviate significantly from the classical result. 
We provide definite evolution equations that can be used to explore potentially important small-scale correc¬ 
tions, such as shock formation, which has been proposed as a mechanism for both burst and quiescent emission 
from magnetars. 


Introduction — From the earliest days of the quantum the¬ 
ory of light, before even the development of quantum elec¬ 
trodynamics (QED) proper, it was recognized that quantum 
effects should become important for electromagnetic field 
strengths of order m 2 /he, where m and e are the mass 
and charge of the electron [1, 2], New effective photon- 
photon interactions emerge, mediated by electron loops, lead¬ 
ing to phenomena such as vacuum birefringence and light-by- 
light scattering. Most dramatically, critical-strength electric 
fields create electron-positron pairs out of the vacuum, a non- 
perturbative effect [3]. The most promising route to reach¬ 
ing these field strengths in the laboratory is the use of high- 
intensity lasers [4]. While some of the effects may be observ¬ 
able in upcoming facilities, the actual field strengths will still 
be subcritic al. 

Fortunately, nature provides us with another avenue to 
investigate strong-field QED: a class of astrophysical ob¬ 
jects known as magnetars. Magnetars are pulsars (rotating, 
magnetized neutron stars) with exceptionally strong surface 
magnetic field strengths. In fact, magnetars can have field 
strengths of up to 10 15 G and maybe higher, which exceed the 
critical field strength, 

o 
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B q = — « 4.4 x 10 13 G, (1) 

he 

by two orders of magnitude! Much work has been devoted to 
understanding the physical processes that take place in such 
magnetic field strengths; see [5-7] for reviews. 

Most of this work is done assuming a vacuum environment, 
whereas in fact magnetars (and pulsars in general) are sur¬ 
rounded by a diffuse plasma. The existence and properties 
of this plasma can be understood from the smallness of the 
dimensionless parameter 


Here B is the magnetic field strength, R is the stellar radius, 
and we have assumed canonical pulsar values B « 10 12 G 
and R ss 10 km. 

This number accounts for the plasma as follows [8, 9], A 
conductor moving with velocity v in a magnetic field B gener¬ 
ates an electric field of order Bv by “unipolar induction”. For 
a rotating magnetized sphere in vacuum this electric field has 
a component along B, and hence can accelerate particles. The 
energy to which the particles can be accelerated over a typ¬ 
ical system size is thus eBvR. Computing v/x shows that 


this energy exceeds the rest mass of the electron by many 
orders of magnitude. (For pulsars a typical surface velocity 
is v ~ 10 -4 .) Thus any stray charges are rapidly acceler¬ 
ated to above the pair-production threshold, and the ensuing 
pair-creation cascade will populate the magnetosphere with 
plasma. 

As charges are generated they will arrange themselves to 
cancel the electric field, driving the Forentz invariant E • B to 
zero. Production ceases as this invariant becomes too small to 
produce the required acceleration. For charge corotating with 
the star the density required to cancel E ■ B is the so-called 
Goldreich-Julian charge density vB/R. This sets a typical 
scale for the plasma mass density, mvB/eR. The ratio of 
the particle mass/energy density to the electromagnetic field 
energy density is then V'x, which is exceedingly tiny, making 
the plasma dynamics completely dominated by the field. 

Assuming classical electrodynamics, such plasmas are de¬ 
scribed by a non-linear theory of the electromagnetic field 
known as force-free electrodynamics (FFE) [8, 10-12], The 
theory follows from the assumption that the electromagnetic 
stress-energy is conserved on its own, leading to the condition 
that the Forentz force density everywhere vanishes. Naively, 
one might expect any classical description to break down at 
or near the critical field strength Bq. The measurement of 
the surface magnetic field strength of a pulsar/magnetar relies 
on the dipole radiation spin-down formula, which has only 
been derived in classical electrodynamics (vacuum [13-15] 
or force-free [16-18]) or in vacuum QED [19]. Thus the very 
evidence for super-critical magnetic fields in nature is sensi¬ 
tive to the question of magnetically dominated QED plasma. 

In this letter we will derive equations of “quantum force- 
free electrodynamics” to describe this plasma. The strategy is 
to integrate out electron loop fluctuations from the QED ac¬ 
tion, following the basic approach established by Euler and 
Heisenberg (EH) in 1935 [2]. However, the EH calculation 
is done assuming no fermion in- or out-states, allowing the 
electron to be integrated out entirely in the effective action, 
whereas we wish to consider plasma. We therefore proceed 
in two steps. First, we consider a collisionless multiparticle 
system and use effective field theory (EFT) arguments to es¬ 
tablish the size and form of the modifications due to QED. We 
then show that the modifications that survive in the magneti¬ 
cally dominated limit follow from the EH Lagrangian, with¬ 
out requiring a new QED calculation. We thereby write down 
definite equations for magnetically dominated QED plasma. 
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Like their classical counterpart, the quantum force-free 
equations imply E • B = 0. This makes the EH Lagrangian 
real and pole-free (no Schwinger pair production), which im¬ 
plies that the effective theory will adequately describe physics 
at and above the critical field strength. The classical force-free 
equations are supplemented by small corrections, and there is 
no drama at the critical scale for magnetars. In particular, we 
expect the classical spindown formula to be corrected only by 
a factor of at most order 10” 3 , supporting the self-consistency 
of the argument for strong magnetic fields in magnetars. 

Small corrections can still have important effects if they 
give rise to qualitatively new features. One new feature here 
is the non-linearity of the vacuum sector of the theory, which 
gives rise to shock formation (e.g., [20]). Heyl and Hernquist 
have argued that such shocks also form in the plasma context 
[21, 22] and that the energy in the shocks and would even¬ 
tually be dissipated into electron-positron pairs, providing a 
mechanism for both burst [23] and quiescent [24, 25] mag- 
netar emission. Their results are based primarily on an anal¬ 
ysis of the characteristics of the one-dimensional, linearized 
problem. Our new, non-linear evolution equations for the full 
three-dimensional description of the quantum plasma provde 
a foundation for exploration of the detailed operation of this 
and other QED effects in realistic magnetosphere models. 

It is noteworthy that a relativistic QED plasma can be given 
such a simple description in the strong-field limit. We show 
how this simplification occurs for any covariant Lagrangian 
assumed to have a standard matter coupling and to have con¬ 
served stress-energy without matter contributions. We discuss 
how every Lorentz-invariant action gives rise to a “force-free” 
dynamics in addition to its ordinary least-action dynamics. 

We use Heavisde-Lorentz units with the speed of light 
c = 1, but leave h ^ 1. Our metric signature is (—, +, +, +). 
We restrict to flat spacetime, but all results generalize straight¬ 
forwardly to an arbitrary curved spacetime. 

Particle EFT — We begin by considering a system of N 
point particles coupled to a gauge field, in the spirit of the 
analyses carried out in Refs. [26, 27]. We will assume that 
these particles do not collide, which is justified by the low 
density of the magnetar plasma. This means that the EFT 
action will be a functional the gauge field A t , together with 
N worldlines zf, with all terms local, and gauge, Lorentz, 
and worldline-reparameterization invariant. If supplemented 
with worldline spinors coupling to the gauge field, this for¬ 
mulation is fully equivalent to ordinary QED [28, 29], except 
that the restriction to fixed N limits its validity to processes 
where electron-position pairs cannot be created. In the deriva¬ 
tive expansion, all effects of the spinor are parametrized by 
their contributions to multipole couplings of the particle. 

As in the derivation of the standard EH Lagrangian, the 
power-counting parameter is the average momentum transfer 
over the particle mass, so that derivatives are suppressed by 
the momentum of the Fourier transform of the field strength 
over the electron mass. In magnetar applications this is the 
electron Compton wavelength over the neutron star radius, 
justifying keeping only to leading non-trivial order in the 
derivative expansion. In practice this means that only terms 
made of A fl , /•],„, and the 4-velocity u 1 ' may appear, and we 


integrate along the worldlines with respect to proper time t, . 
The complete collection of terms is 1 


Seff = 


£em [Fnu\d x 


A* 

N 

E 

t=i 


[ - qi u^A^) dn. (3) 

J Zi 


Here is a constant interpreted as the particle charge. If rri, 
were likewise constant, it would be the standard notion of par¬ 
ticle mass. However, here we say only that rrii is a function 
with dimensions of mass that is constructed covariantly and 
gauge-invariantly from F t , v and u 1 -'. Finally Lem is an arbi¬ 
trary function of F l: v only, and hence can depend only on the 
invariants 


I = F pv F^, K = F^F^, F^ = \e^ pa FPF (4) 

In (3 + 1) langauge these invariants are I = 2(1 F — E 2 ) and 
K = 2E • B. The “mass” rn, can depend on these invariants 
as well as the third invariant F^FP pu a u^. 

Varying with respect to the gauge field A 1 ' gives 

+ = JA1 (5) 


with 


N 


= E / ® u ' 1 + 2V - 


»=i ■ 


drrii 

~dFuv 


S(x — Zi(r))dTi. (6) 


Varying with respect to z M gives for each particle i. 


■^T \ ™,iU f + ty 1 " + <<) 


dnii 


— Qi Fu 


drrii 


dF c 


Oifi 


-V^Fa/? (7) 


As with ordinary electromagnetism coupled to point particles, 
Eqs. (5)-(7) mix (5-functions with nonlinearities and hence 
must be regarded as formal. The origin is that the EFT should 
not describe physics on the scale of the particle size. Re¬ 
stricting to the appropriate larger scales, each particle makes a 
small perturbation to the long-wavelength field produced by 
the other particles, and the delta functions appear as point 
sources for the linearized version of Eq. (5). Equivalently, one 
can think of the force law (7) as holding for the field produced 
by all other charges except the one under consideration. The 
neglected self-force effects are small in the magnetar applica¬ 
tion, but may be important for laser plasmas [4]. 


1 Dipole couplings also appear at this order in the derivative expansion, and 
technically should be included in Eq. (3). However, after matching to QED 
and taking the magnetically dominated limit these terms are the same order 
as those we drop. For simplicitly we do not write them down at this stage. 

2 Since depends on the worldline, one should instead vary with respect to 

z? (A) for some fixed parameter A. 
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QED Matching and the Magnetic Limit — The unknowns 
qi, rrii, and are to be matched with a UV the¬ 

ory valid on small scales, in this case QED. The charges qi 
match trivially with ±e (with + for positrons and — for elec¬ 
trons) based on the form of the interaction term in the QED 
Lagrangian. To determine Lem we note that this term does 
not depend on u 1 ' and hence can be considered in the spe¬ 
cial case N = 0 where there are no particles in the EFT. 
On scales smaller than the typical variation R of the field the 
field strength F^ v is approximately constant. Thus we want 
to match to QED in the presence of a uniform external field. 
Since there are no fermion degress of freedom in the EFT we 
must integrate these out of QED. This is the classic Euler- 
Heisenberg calculation, and hence 

£em = £eh, (8) 


where the expression for Leh may be found in standard texts 
(e.g., Eq. (1.2) of Ref. [30] orCh. 33 of Ref. [31]). 

The functional form of the “mass” m* [F^ v , u^] will be pre¬ 
sented in a future paper. 3 For the present paper we will need 
only the fact that m, is a real and slowly-growing function of 
field strength in the magnetic case I > 0, K = 0 of relevance 
here, which follows from results already known in the liter¬ 
ature. In particular, the electron ground state energy shift in 
a strong magnetic field [34, 35] is order am in all interesting 
magnetic field strengths. There are also corrections to rrii pro¬ 
portional to the velocity it M , but these follow from matching 
to the same calculation in the UV theory, namely, the renor¬ 
malization of the exact electron propagator in a background 
field, and so must have the same analytic structure. 

The (related) facts that the magnetar field is magnetic (I > 
0) and the effective mass is order m turn out to allow us to 
drop all terms involving m, for that application. For exam¬ 
ple, in Eq. (6) the spacetime derivative counts as a power of 
inverse length 1 /R, while the field strength derivative counts 
as a power of 1 /B. Thus dividing the second term by the first, 
the ratio m/eBR is the parameter y defined earlier in Eq. (2), 
and the current is given to an excellent approximation as 



qiU^S(x 


Zi{r))dTi. 


(9) 


The mass terms in the force law (7) are also formally su- 
pressed in this way, but dropping them requires more care, 
and relies crucially on the field being magnetic, rather than 
electric or null. The subtlety is that the LHS of Eq. (7) con¬ 
tains terms of higher differential order, which therefore can 
have their scale set by the dominant (Lorentz) force together 
with initial conditions. In the magnetic case this scale is the 
cyclotron motion about the field line, and dropping m from 
the Lorentz force law averages over the gyrations [36, 37]. 
(In strong fields the cyclotron motion is quantized, so energy 


3 This term gives rise to strong-held corrections to the Lorentz force law, of 
the variety considered in Ref. [32, 33]. 


eigenstates don’t show any actual gyration.) Setting rrii = 0 
in Eq. (7), we arrive at 

= 0 . ( 10 ) 

The content of this statement is that the particles are stuck to 
magnetic field lines, moving freely along them in the guiding 
center approximation. For further discussion of the relativistic 
meaning of this statement, see Sec. 3.2.3 of [12], 

The relationship between the small parameter y and the 
EFT (3) is clarified by rewriting the former as 

m Bn h 

y =- =-B- - (ll) 

A eBR B mR v ’ 

The second factor ( h/m)/R is the power-counting parameter 
of the derivative expansion. Thus y being small in large mag¬ 
netic fields automatically implies the worldline EFT is valid. 
On the other hand, the EFT can still be valid when y is large, 
such as if the fields are nearly null (B = s/I/2 -C Bq) but 
still strong, F IIV ~ Bq. The derivation above emphasizes the 
separate physical origin of the quantum corrections from the 
mass-independence of the particle equations of motion, high¬ 
lights the role of the magnetic assumption in the latter, and 
avoids subtleties with the orders of limits involved in keeping 
to finite to derive the EH terms, while setting to = 0 in the 
particle worldlines. 

Quantum FFE — Dotting both sides of Eq. (9) with 
and using Eq. (10) gives the force-free condition. 


F )lv .r = 0 . ( 12 ) 

This is the statement that the Lorentz force density vanishes 
at every point in the plasma. Using (5) to eliminate J M in 
favor of we obtain a closed system of equations for the 
electromagnetic field, without reference to the charges. (The 
charges move on field lines to provide the current, but their 
mass does not affect the field dynamics.) 

The existence (12) of a zero-eigenvector F' for F vv implies 
that the second invariant vanishes, K = 0. Such fields are 
called degenerate and have a beautiful mathematical structure 
[10, 12, 38]. They define a foliation of spacetime by two- 
surfaces, which in the magnetic case I > 0 are timelike and 
represent worldsheets of magnetic field lines, or field sheets. 
They can always be expressed as the wedge of two one-forms 
F/j,v = O’hj.b^}, which themselves are exact, a M = and 

bfj. = V^ 2 for two scalar “Euler potentials” 0i and 02- 
In light of the degeneracy of F IIV we may now simplify 
the force-free equations. Substituting Eq. (5) into Eq. (12) 
and using F^ v = it follows that the term involving F 

always vanishes, and we obtain 

F/ivS/p {f(I)F pv } = 0, /=- . (13) 

(JI K =0 

Using the explicit form of the EH Lagrangian, we have 

/ = ! + <? + 0(a 2 ), (14) 


with 



( 15 ) 
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Here we introduce the dimensionless magnetic field strength 
b = B/Bq = \jil/‘2 /Bq. The integral in (15) is finite and 
can be computed numerically at a given value of b or expanded 
at strong or weak fields [30]. It is of order 10~ 3 at magnetar 
field strengths. This expression includes all terms involving 
one fermion loop. It is corrected at 0{a 2 ) by higher loops, 
but is non-perturbative in the field strength b, holding at arbi¬ 
trarily 4 strong fields. Eqs. (13), (14) and (15) are the theory of 
quantum force-free electrodynamics in the one-loop approxi¬ 
mation. 

Evolution form — While the elegant form (13) is conve¬ 
nient for general manipulations, for numerical solution one 
must express the theory in the form of evolution and constraint 
equations. We take the electric and magnetic fields E and B 
as our fundamental variables. We begin by defining a.fictitious 
current J 1 ' by 

( 16 ) 

in terms of which the equations of motion (13) still take the 
form F t , v ,]>• = 0. In this form it is clear that the inho¬ 
mogeneous Maxwell equations are modified by E —> /E, 

I 


B —> f B, p —> p and J — > J, 

V • (/E) = p (17) 

Vx(/B)=j + ft(/E). (18) 

The homogeneous Maxwell equations V [ /( F pcr \ = 0 are un¬ 
modified, 

V • B = 0 (19) 

V x E = -d t B. (20) 

The force-free condition F. JV J 11 0 is equivalent to 

E • J = 0, pE + J x B = 0. (21) 

The fact that K = 0 for any solution is equivalent to E B = 0, 
whose time derivative states that 

E • B = 0 <9 t E ■ B = E ■ (V x E), (22) 

where (20) has been used. 


Eq. (20) is the evolution equation for B. We determine the evolution equation for E by projecting <9 f E onto the basis 
{E, B, E x B}, using Eqs. (17)-(22), and eliminating p and J. The result is 




E f E • V x (/B) + 4f'E 2 B • V x E 


/ - 4/'£ 2 


Jr {E ■ v x E} + {ExB-Vx (/B) - E 2 V ■ (/E)} , 


E 2 B 2 / 


(23) 


where a prime represents an /-derivative. If we linearize with respect to g then we find 


<9 t E = V x B + — (E - VxE — B-VxB) — V-E 
B- 


E x B 

B 2 


dg E x B 
dl B 2 


E-V/+— E 

dl 


ExB-V/ 

E 2 


4(E - V xB + B- V xE) 


(24) 


where we remind the reader that / = 2 (B 2 — E 2 ). The first line gives the evolution equation for standard Maxwellian FFE with 
the last two terms comprising (minus) the Maxwellian current J. We may view the remainder of Eq. (24) providing effective 
corrections to the current. The first term of the second line of Eq. (24) corrects the “advection" term V ■ E (E x 14 / B 2 ) = pv 
(with p the charge density and v the drift velocity) for the effective charge density (17). The second term introduces an effective 
current directed along the electric field, a qualitatively new feature. (The real current is orthogonal to the electric field, as 
required by the force-free equation ( 12).) 

The evolution system consists of evolution equations <9,B = —V x E and either Eq. (23) or (24) for 9 t E, together with 
constraints V • B = 0 and E • B = 0. By construction, the evolution equations preserve the constraints. We have not determined 
whether this formulation is mathematically well-posed. If g = 0 the sytem is equivalent to ordinary FFE, which can be made 
hyperbolic provided I > 0 [39, 40]. 


Previous Work — Heyl and Hernquist postulated the La- 
grangian Leu + OK, where 0 acts as a Lagrange multiplier. 
The equations of motion are then K = 0 as well as Eq. (5) 
together with the expression -I 1 ' = 0,JJ F I “ J for the current. 5 


4 Since £em grows only like a log b at large b. there is no reason to suspect 
a breakdown of perturbation theory until the fantastical value of b ~ e 137 , 
corresponding to field strengths of ~ 10 60 G. 

5 Thompson and Blaes [41] have also argued directly for this form of the 


They perturbed from a uniform field and eliminated the La¬ 
grange multiplier in this (1 + 1)D linearized theory. The La¬ 
grange multiplier may be eliminated in general by recalling 
that K = 0 implies F jn/ = for some (non-unique) pair 

of one-forms. Dotting with the above expression for the cur¬ 
rent gives the force-free condition F tiy J u = 0. Conversely, 


current based on bozonization of the fermion field on each field line. 
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one may show that the force-free condition implies that form 
for some scalar 9. Thus the Lagrangian gives rise to the full 
(3 + 1), non-linear theory we consider. 

Generalized force-free fields — We conclude with a dis¬ 
cussion of force-free theories generally. Consider a La¬ 
grangian L assumed only to be a covariant, gauge-invariant 
functional of the metric and gauge field, 

S = J L [g^,A ll \^gd 4x . (25) 

Instead of imposing the equations of motion we define their 
failure to be satisfied to be the stress-energy and charge- 
current, 

rp _ _2_i^ v , _ _ 

y/=g8g^' y/=g8A„ ( } 

The first is the standard definition of stress-energy, while the 
second is equivalent to adding a coupling term to the 

action (25). The covariance and gauge-invariance of Lem im¬ 
ply that these quantities satisfy (e.g., [42]) 

= 0, V"T /t „ = -F^r. (27) 

No field equations are imposed here; these are simply iden¬ 
tities that hold for arbitrary vector fields A 1 11 '. We give the 
theory content by demanding conservation of stress-energy, 
t „ = 0, which immediately implies the force-free con¬ 
dition F IJtl ,J v = 0. Allowing non-zero current corre¬ 
sponds to allowing unspecified coupling to matter, but de¬ 


manding conservation of l ' llv consists of neglecting the stress- 
energy of the matter. This is appropriate when there is much 
more energy in the field than in the matter. Thus we pro¬ 
vide a simple derivation that the strong-field limit is force- 
free for an arbitrary classical electrodynamics (such as Born- 
Infeld). It would be an efficient derivation of the force-free 
equations for QED plasma if the coupling A^ and assump¬ 
tion = 0 could be justified directly from QED in the 

appropriate limit. 

This type of calculation need not be restricted to electrody¬ 
namic theories, and indeed an identity of the form = 

Y[ip,8S/5ip] exists for every covariant Lagrangian of some 
collection of fields t/> [42-44]. By demanding = 0 

one can define “force free” conditions Y = 0 for any such La¬ 
grangian. Thus, besides the usual field equations 5S/8f> = 0, 
each covariant Lagrangian gives rise naturally to a second set 
of equations Y[ip, 5 S/Si/j\ = 0. 
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